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I. INTRODUCTION 



The bound and scattering states of the s- and /-waves for any interaction system have 
raised a great interest in non-relativistic as well as in relativistic quantum mechanics [1-3]. 
The exact solution of the wave equation is very important since the wavefunction contains 
all the necessary information regarding the quantum system under consideration. A number 
of methods have been used to solve the wave equations exactly or quasi-exactly for non-zero 
angular momentum quantum number (I ^ 0) by means of a given potential. The bound 
state eigenvalues were solved numerically [4,5] and quasi-analytically using variational [4,6], 
perturbation [7], shifted 1/N expansion [8,9], NU [10,11], SUSYQM [12-14] and AIM [15] 
methods. 

The Hulthen potential [10,12,13,15,16] is one of the important short-range potentials in 
physics and it has been applied to a number of areas such as nuclear and particle physics 
[17], atomic physics [18,19], molecular physics [20,21] and chemical physics [22]. Therefore, 
it would be interesting and important to solve the relativistic equation for this potential for 
/ 7^ case since it has been extensively used to describe the bound and continuum states of 
the interaction systems. Recently, the exact solutions for the bound and scattering states 
of the s-wave Schrodinger [16,23], Klein-Gordon [1-3] and Dirac equation [24,25] with the 
scalar and vector Hulthen potentials are investigated. 

Relativistic effects with the scalar plus vector Hulthen- type potential [1,2] in three- and 
D dimensions and harmonic oscillator-type potential [26,27] have been also discussed in the 
literature. The bound-states of the Dirac and Klein-Gordon equations with the Coulomb-like 
scalar plus vector potentials have been studied in arbitrary dimension [28-32]. Furthermore, 
the exact results for the scattering states of the Klein-Gordon equation with Coulomb-like 
scalar plus vector potentials have been investigated in an arbitrary dimension [33]. This 
equation has been exactly solved for a larger class of linear, exponential and linear plus 
Coulomb potentials to determine the bound state energy spectrum using two semiclassical 
methods with the following relationship between the scalar and vector potentials: V(r) = 
V + (3S(r), S(r) > V(r) where V and (3 being arbitrary constants [34]. In particular, 
inserting the constants Vq = and (5 — ±1 provides the equal scalar and vector potential 
case V(r) = ±S(r). 

Also, the position-dependent mass solutions of the nonrelativistic and relativistic systems 
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have received much attention recently. Many authors have used different methods to study 
the partially exactly solvable and exactly solvable Schrodinger, Klein-Gordon and Dirac 
equations in the presence of variable mass having a suitable mass distributions function 
in ID, 3-D and/or any dimension D cases for different potentials, such as the exponential- 
type potentials [35], the Coulomb potential [36], the Lorentz scalar interactions [37], the 
hyperbolic- type potentials [38], the Morse potential [39], the Poschl- Teller potential [40], the 
Coulomb and harmonic potentials [41], the modified Kratzer-type, rotationally corrected 
Morse potentials [42], Mie-type and pseudoharmonic potentials [43]. Recently, the point 
canonical transformation (PCT) has also been employed to solve the .D- dimensional position- 
dependent effective mass Schrodinger equation for some molecular potentials to get the exact 
bound state solutions including the energy spectrum and corresponding wave functions [41- 
43]. 

A new method to obtain the exactly solvable PT-symmetric potential potentials within 
the framework of the variable mass Dirac equation with the vector potential coupling scheme 
in (1 + 1) dimensions [38]. Three PT-symmetric potentials are produced which are PT- 
symmetric harmonic oscillator-like potential, PT-symmetric of linear plus inversely linear 
potential and PT-symmetric kink-like potential. The SUSYQM formalism and function 
analysis method are use to obtain the real energy levels and corresponding spinor components 
for the bound states. Further, the position-dependent effective mass Dirac equation with 
the PT-symmetric hyperbolic cosecant potential can be mapped into the Schrodinger-like 
equation with the exactly solvable modified Poschl- Teller potential [38] . The real relativistic 
energy levels and corresponding spinor wavefunctions for the bound states have been given 
in a closed form. 

The Nikiforov-Uvarov (NU) method [44] and other methods have also been used to solve 
the .D- dimensional Schrodinger equation [45] and relativistic D-dimensional Klein-Gordon 
[46] , Dirac [47] and spinless Salpeter equations [48] . 

In strong coupling cases, it is crucial to understand relativistic effects on a moving par- 
ticle in a potential field. In a non-relativistic case, Schrodinger equation with the Hulthen 
potential [10,12,13,15] was solved using the usual existing approximation, 4y ps a 2 ^Sl^ 
for the centrifugal potential which was found to be consistent with the results of other 
methods [4,8,13,15]. Unfortunately, this approximation is valid only for small values of the 
screening parameter a, but the agreement becomes poor in the high-screening region [10, 15]. 
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Hence, it is of sufficient need to improve the analytical results for the Schrodinger equation 
with the Hulthen potential by means of a new approximation scheme. Recently, Haouat 
and Chetouani [49] have solved the Klein-Gordon and Dirac equations in the presence of the 
Hulthen potential, where the energy spectrum and the scattering wavefunctions are obtained 
for spin-0 and spin-| particles, using a more general approximation scheme, ^ ~ <a 2 ^ e J_^ r ^ 
where 7 is a dimensionless parameter (7 = 0, 1 and 2) for the centrifugal potential. They 
found that the good approximation, however, when the screening parameter a and the di- 
mensionless parameter 7 are taken as a = 0.1 and 7 =1, respectively, which is simply the 
case of the usual approximation [10,12,13,15]. Also, Jia and collaborators [50] have recently 
proposed an alternative approximation scheme, ^ ~ a 2 (^ ea r_ 1 + wnere u is a di- 

mensionless parameter (uj = 1.030), for the centrifugal potential to solve the Schrodinger 
equation with the Hulthen potential. Taking uj — 1, their approximation can be reduced 
into the usual approximation [10,12,13,15]. However, the accuracy of their numerical results 
[50] is found to be in poor agreement with the other numerical methods like integration 
and variational methods [4,5]. In order to improve the accuracy of the used approxima- 
tion, we propose and apply an alternative shifted approximation scheme to approximate the 
centrifugal term given by [51,52] 
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where Co is a shifting dimensionless parameter. The approximation scheme (1) emerged as 
a quite successful formalism to study the Schrodinger equation with the Manning-Rosen, 
hyperbolic and Hulthen potentials in calculating the energy eigenvalues within the frame- 
work of the NU method [51-53]. The accuracy of the results are significantly improved 
over all other existing literature approximation schemes and analytical methods [13,15,50]. 
With extremely high accuracy, we have obtained the numerical energy eigenvalues as with 
those obtained by the numerical integration [4,5,53], variational [4] methods and also by a 
MATHEMATICA package programmed by Lucha and Schoberl [54]. 

The purpose of this work is to employ the approximation scheme given in (1) to solve 
the position-dependent mass radial Klein-Gordon equation with any orbital angular quan- 
tum number I for the scalar and vector Hulthen potentials in D-dimensions. This offers a 
simple, accurate and efficient scheme for the exponential-type potential models in quantum 



4 



mechanics. 

Our paper is organized as follows. In section 2, we review the NU method. In section 3, 
we present a brief a derivation to find the shifting parameter c . Then, the analytical solution 
of the position-dependent mass Klein-Gordon equation with the scalar and vector Hulthen 
potentials is obtained for any /-state by means of the N-U method. Section 4 contains the 
summary and conclusions. 



II. NU 



The NU method is breifly outlined here and the details can be found in [44] . This method 
is proposed to solve the second-order differential equation of the hypergeometric type: 

r n (z) + ^<{z) + £j^n(z)=0, (2) 

where a(z) and a(z) are polynomials, at most, of second-degree, and r(s) is a first-degree 
polynomial. In order to find a particular solution for Eq. (2), let us decompose the wave- 
function ip n (z) as follows: 

i/> n (z) = (j) n (z)y n (z). (3) 
We can reduce Eq. (2) into the form 

a{z)yl{z) + r(z)y' n (z) + Xy n (z) = 0, (4) 

with 

t(z) = t(z) +2tt(z), t'{z) <0, (5) 
where r'(z) = ^j^- is the derivative. Also, A is a constant given in the form 

\ = \ n = -nr\z)-^n(n-l)a"(z), n = 0, 1, 2, • • • , (6) 

where 

\ = k + 7r'{z). (7) 
The y n (z) can be written in terms of the Rodrigues relation 

where B n is the normalization constant and the weight function p(z) satisfies 

a(z)p'(z) + (a'(z)-r(z))p(z) = 0. (9) 
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The other wavefunction in the solution is defined by 

a(z)<f)'(z) -ir(z)<f>(z) = 0. (10) 
Further, to find the weight function in Eq. (8) we need to obtain the following polynomial: 

n(z) = \ \a\z) - r(z)] ± ji [a\z) - r(z)] 2 - a(z) + ka(z)^ . (11) 

The expression under the square root sign in Eq. (11) can be arranged as the square of a 
polynomial. This is possible only if its discriminant is zero. In this regard, an equation for 
k is being obtained. After solving such an equation, the determined values of k are included 
in the NU method. 



III. BOUND-STATE SOLUTIONS 



A. An Impoved Shifted Approximation Scheme 



The approximation is based on the expansion of the centrifugal term in a series of expo- 
nentials depending on the intermolecular distance r. Therefore, instead of using the usual 
existing approximation in literature, let us, instead, take the following exponential-type 
potential to approximate the centrifugal potential, 



1 



a 2 [co + v(r) + v 2 (r)] , v(r) 



-1' 
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(e ar - 1) 



(12) 



In the low-screening region, 0.4 < ar < 1.2 [15] (i.e., small screening parameter a), Eq. (12) 
is a very well approximation to the centrifugal potential and the Schrodinger equation for 
such an approximation can be easily solved analytically. In Fig. 1, we give a plot of the 
variation of the centrifugal potential and its approximation given in Eq. (12) versus ar. It 
shows that the approximation (12) and 1/r 2 are similar and coincide in both high-screening 
as well as in the low-screening regions. 

Changing the r coordinate to x by shifting the parameters as x = (r — r )/r to avoid 
singularities [55], we obtains 

1 1 



\(l + x)~ 2 



a 



c + 



+ 
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7 = tt^o, 



(13) 
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and expanding Eq. (13) around r = r (x — 0), we obtain the following expansion: 



l-2x + 0{x 2 ) = 7 2 c + 



1 



+ 



1 
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eT - 1 ( e 7 - \y ( e 7 - 1) 



a; + 0(a; 2 ), 



(14) 



and consequently 
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c + 
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3 2 



+ 



2. 



,eT-l ' (eT-1) 2 (eT-1) 3 
By solving Eqs. (14) and (15) for the dimensionless parameter c , we obtain 

1 1 1 



c 



r 



eT - 1 (eT - l) 2 



0.0823058167837972, 



(15) 



(16) 



where e = 2.718281828459045 is the base of the natural logarithms and the parameter 
7 = 0.4990429999. 

Therefore, the centrifugal potential takes the form 

2 



lima 2 

a^O 
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r 



+ 



eT - 1 (eT - l) 2 1 - e 



+ 



1-e" 



(17) 



Let us remark at the end of this analysis that the approximation used in many papers in 
literature [10,12,13,15] is a special case of Eq. (12) if c is set to zero. 



B. A Quasi-Exactly Energy Eigenvalues and Eigenfunctions 

The D-dimensional time-independent radial position-dependent mass Klein-Gordon equa- 
tion with scalar and vector potentials S(r) and V(r), respectively, r = |r| , and position- 
dependent mass m(r) describing a spin-zero particle takes the general form [3,46] 

V^^M + ^2" {[Em ~ V(r)f - [m(r)c 2 + S(r)] 2 } V&JW = 0, 

V ^ = E W = Ri{r)Yti D -^ 02,-, On-,), (18) 

where E n i denotes the Klein-Gordon energy and V 2 D denotes the D-dimensional Laplacian. 
Further, x is a D-dimensional position vector. Let us decompose the radial wavefunction 
Ri(r) as follows: 

Rt(r)=r-V>-Wg(r), (19) 
7 



we, then, reduce Eq. (18) into the following D-dimensional radial position-dependent effec- 
tive mass Schrodinger-like equation 



d 2 g(r) , 1 



dr 2 h 2 c 



2 (D + 21 - 1)(D + 21 - 3)h 2 c 2 



[E nl -V(r)] 2 -[m(r)c 2 + S(r)Y- 



4r 2 



Further, taking the vector and scalar potentials as the Hulthen potentials 



V(r) 



V e 



1 - e- ar 

and choosing the following mass function 



, S(r) 



S e~ 



1-e- 



a = r o > 



g(r) = 0. 
(20) 



(21) 



we can rewrite Eq. (20) as 



9"(r) + 



m(r) =m + 

1 — e 



1 / 2 [m c 2 (S-q - m lC 2 ) + £„jVb] e" 



(22) 



h 2 c 2 { i. - ( 

2\2] „-2ar ft 2 c 2 a 2 , 



Vtf - (^o - m lC 2 r e~ 2OT - ^L[D + 2l-l)(D + 2l- 3)e 



(1 - e~ ar Y 



g(r) 



h 2 c 2 



(m c 2 ) 2 - E 2 nl + AE t g(r), g(0) = 0, 



(23) 



with the shift energy AE t = h 2 c 2 a 2 (D + 21 - l)(D + 21 - 3)c /4. On account of the wave 
function g(r) satisfying the standard bound-state condition (real values), i.e., g(r — > oo) — > 
0. If we rewrite Eq. (23) by using a new variable of the form z = e~ ar (r e [0, oo), z G [1, 0]), 
we get 



d 2 g(z) 1-z dg(z) \ 



dz 2 ' z(l-z) dz ' \z(l-z)} 2 
x {-4 + (A - A - 7 + 2e*,)s - (ft + ft + ft - ft + e 2 ,) s 2 } g{z) = 0, 

where the following definitions of parameters 
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\l (™oc 2 ) 2 - E 2 nl + AE t _ 2 ( moC 2s + EnlVo ) 



Q 



, A 



-, A 



o2 _ T/2 



mic 2 (mic 2 — 25o) 2momic 4 
A = 7^ » A = 7v? > 7 



g 2 ' " g 2 

(D + 2Z- l)(D + 2Z-3) 



g 2 



g 2 



, g = ftca, (25) 
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are used. For bound-state solutions, we require that V < (So — mic 2 ) and E n i < 
a/ (m c 2 ) 2 + AEi. In order to solve Eq. (24) by means of the N-U method, we should 
compare it with Eq. (2). The following values for parameters are found 



t(z) = 1-z, a(z)=z-z 2 , a(z) = -e 2 nl +{(3 l -(3 A - 1 +2e 2 nl )s-{(3 l +(3 2 +^-(3 A +e 2 nl )s 2 . (26) 

If we insert these values of parameters into Eq. (11), with cr\z) = 1 — 2z, the following 
linear function is obtained 



= -f ±2 \fc + 4 (A + & + ft - ft + zli - k)] z 2 + [4(fc - ft + ft + 7 - 2^,)] z + 44 

(27) 

The determinant of the square root must be set equal to zero, that is, A = (k — ft + ft + 
7 - 2e 2 nl ) 2 - e 2 nl [1 + 4(ft + ft + ft - ft + £„, - k)] = 0. Thus, the constant k found to be 



k = ft - ft - 7 ± £„/ v / l + 4(ft + ft + 7 )- 

In this regard, we can find four possible functions for ir(z) as 

z , I e n i =f [e„i - | VI + 461 z for fcj = d + e nX ^\ + 46, 
7r(s) = --±< 

2 [ =F + |Vl + 46J 2 for fc 2 = d - fWl + 46. 

where 6 = ft + ft + 7 and d = ft — ft — 7. Thus, taking the following values 



fc = ft - ft - 7 - <Wl + 4(ft + ft + 7), 



(28) 



(29) 



(30) 



and 



they give 



n(z) = -- +e ni 



Snl + ^1 + 4(^2 + ^3 + 7) 



Hz) = 1 + 2e. 



ni 



l+£n« + ^V /l + 4 (ft + ft + 7) 



r s 



1 + £n« + ^1+4(^2+^3 + 7) 



< 0. 



(31) 



(32) 



Eqs. (30)- (32) together with the assignments given in Eq. (26), the following expressions 
for A are obtained 



A n = A = n 2 + 



1 + 2e nl + v 7 ! + 4(ft + ft + 7 )1 n, (n = 0, 1, 2, • • • ), 



A = ft - ft - 7 - -(1 + 2e nl ) [l + v/l + 4(ft + ft + 7 ) 



(33) 
(34) 
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where n is the radial quantum number. By defining 



S = \ (l + ^1+4(^2 + Js + 7)) , (35) 

where f3 2 + 13% = 5 2 — 5 — ^ . With the aid of Eq. (35), we can easily obtain the energy 
eigenvalue equation of the Hulthen potential by solving Eqs. (33) and (34): 

( D)_ ((3 1 -(3 4 - 1 -n 2 )-(2n+l)5 



"nl 



2(n + 5) 

4 - & - n 2 - (2n + 1)5} -(D + 21- 1)(D + 2/ - 3) 



8(n + 5) 



m c 2 S + E±M 



+ S 2 



V 2 



n + 5 



(n = 0,l,2,.-.), 



(36) 



2Q 2 (n + 5) 2 
where S = So — mic 2 is the modified scalar potential. Solving the last equation for the 
energy eigenvalues we obtain 



nl 2 



4S S + 2m c' 



W 2 + K, 



nl 



_|_ K nl 

~2~ 



\ 



4S S + 2m c< 



-i 2 



4F 2 + «; 



nl 



(2m c 2 ) 2 + h 2 c 2 a 2 (D + 2l-l)(D + 2l- 3)c 



4V 2 + K. 



= ftcct (2ra + 1) + W4 (S6 2 - \/ 2 ) + (ftca) 2 (£> + 21 - 2) 2 , 



(37) 



0,1,2, 



where n = 0, 1, 2, • • • and / 
quantum numbers, respectively, and 



signify the usual radial and angular momentum 



(hca) 2 (D + 21- 2) 2 + 4S 2 > W 2 , 4£ > 



1 - 



4S [So + 2m C 
4K 2 + «2 



-i 2 



(38) 



are constraints over the strength of the potential coupling parameters. In the above equation, 
let us remark that it is not difficult to conclude that all bound-states appear in pairs, two 
energy solutions are valid for the particle E p = E^ and the second one corresponds to the 
anti-particle energy E a = E~ t in the Hulthen field. When we take the scalar and vector 
potentials as S = (i.e., S = mic 2 ) and V ^ 0, the energy equation (37) becomes 



E nl 



V 1± Kni l(2m c 2 ) 2 + h 2 c 2 a 2 (D + 2l-l)(D + 2l-3)co 1 
2 2 



W 2 + K 



nl 



+ 
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(4m c 2 ) 2 + Ah z c 2 a\D + 21 - 1){D + 21 - 3)c > AV 2 + <„ 



Knl = hca(2n + 1) + \J(hcaY (D + 21 - 2) 2 - 4\/ 2 , D > 1, 
with the following constraints on the coupling parameter of the vector potential: 

(hcaf (D + 21 — 2f > 4V 2 , 



(39) 



(40) 



must be fulfilled for real eigenvalues. 

Therefore, having solved the D-dimensional position-dependent mass Klein-Gordon equa- 
tion for scalar and vector usual Hulthen potentials, we should make some useful remarks. 

(i) For s-wave (/ = 0), the exact energy eigenvalues of the ID Klein-Gordon equation 
becomes 



n ~ Y 



AS S + 2m c 2 



4 V 2 + kI 



i for. 




ASv (S + 2m c 2 



4V 2 + kI 



n n = hca (2n + 1) + J (hca) 2 + 4(S 2 - V 2 ), 



(41) 



In order that at least one level might exist, it is necessary that the inequalities 

,,,„ . I / IS, (,S' (1 + 2,„ (1 r 2 ) \ " 



h z c z a z + 4S 2 > 4K 2 , 



16mQC 4 
4K, 2 + nl 



> 



W 2 + nl 



(42) 



are fulfilled. In the case So = 0, Vq ^ 0, the energy spectrum (in units where H — c — 1): 
Vb 



El 



± 



a (2n + 1) + yja. 



2 - 4K, 2 



ms 



\ 4V 2 + [a (2n + 1) + y/a 2 - AV 2 



1 

16' 



(43) 



with the following constraints on the potential coupling constant: 



16m 2 > W 2 



^/a 2 -AV 2 + {2n + l) 



, a > 2V , 



(44) 



are fulfilled for bound state solutions. We notice that the result given in Eq. (43) is identical 
to Eq. (31) of Ref. [56]. As can be seen from Eq. (43), there are only two lower- lying states 
(n = 0, 1) for a Klein-Gordon particle of rest mass mo = 1 and screening parameter a — 1 
with vector coupling strength Vq < 1/2. As an example, one may calculate the ground state 
energy for the vector coupling strength V = a/2 as 



Vo 
2 



1 ± 



'2m 2 



(45) 
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Further, in the case of pure scalar potential (Vq = 0, Sq = m\c 2 \ the energy spectrum 



Et = ±\l™y - {hCa) + , 4m 2 c 4 > (hcaf (n + I) 2 



(46) 



Since the Klein-Gordon equation is independent of the sign of E n for scalar potentials, the 
wavef unctions become the same for both energy values. If the range parameter a is chosen 
to be a = 1/A C , where A c = h/rriQC = l/m denotes the Compton wavelength of the Klein- 
Gordon particle. It can be seen easily that while So — > mic 2 in ground state (n — 0), all 
energy eigenvalues tend to the value E ss 0.866 m . 

(ii) For D = 3, the mixed scalar and vector Hulthen potentials, the energy eigenvalues 
for I 7^ are given by 



E ± = Y± 

nl 2 



' 4S (s 



+ 2m c 2 



W 2 + kI 



± K 



nl 



AS (S + 2m c 2 ) 



(m c 2 ) 2 + fi 2 c 2 a 2 l(l + l)c 
4^ + ^ 



/v„, = //m(2/, 1) (hcaf (21 • 1 ) 2 • 1 (s 2 V 2 ). 



(47) 



Further, in order that at least one real eigenvalue might exist, it is necessary that the 
inequality 

' 4S [[So + 2m c 2 ) ^ 2 



(ftca) 2 (2/ + l) 2 + 4S 2 > 4^ 2 , 4£ > 



4^ 2 + Ki 



(48) 



must be fulfilled. For the case where S = in the spatial-dependent mass (So = 0, in the 
constant mass case) [46], the energy eigenvalues turn out to be 



± _ Vo l (m c 2 ) 2 + h 2 c 2 a 2 l(l + l)c _ J_ 

^-2 ±rinl ]j W 2 + V 2 nl 16' 

rj nl = hca (2n + 1) + \J (hcaf {21 + l) 2 - W 2 , hca(2l + 1) > 2V , (49) 
with the following constraint over the potential parameters: 

(4m c 2 ) 2 + 16h 2 c 2 a 2 l(l + l)co > 4\/ 2 + hca{2n + 1) + yj {hcaf {21 + l) 2 - 4V 2 . (50) 
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(iii) When D = 3 and I = 0, the centrifugal term { - D+21 1 ^ +21 3 - ) = and consequently the 



approximation term (D+2 ^ 1)( f +2 '~ 3)a2 
ues turn to become 



c ° + 71 



= 0, too. Thus, the energy eigenval- 



v/(m c 2 ) 2 - £± 2 = 



m c 2 S + E±V 



+ s 2 ~v 2 



2hca(n + 5) 



hca 



n + 5 



'-5 



1 + 



(hca) 



(hca) 



V 2 



, (n = 0,l,2,3,..-) 



(51) 



which gives 



E. 



± 



Vo ( AS (S + 2m c 2 ) 

~ 1 1 W^Tl 



(rriQC 2 )' 



\ W 2 + , 2 4 



1 - 



4^0 (S + 2m c< 
4V 2 + <* 



S n = hca (2n + 1) + ^J(hcaf + A(S 2 -V 2 ^, 



(hca) 2 + AS 2 > 4V 2 , (4m c 2 ) 2 > (4V 2 + 



AS (s + 2m c 2 ) 



(52) 



(iv) For equal scalar and vector usual Hulthen potential (i.e., So = Vo), Eq. (36) with the aid 
of Eq. (25) can be reduced to the relativistic energy equation (in the conventional atomic 
units h = c = 1): 



2 (D + 21- 1)(D + 21-3)cq 
m ° + 4T 2 



El 



2r V [m + E R — m 1 \+r (m 1 — 2rn )rn 1 n + 5 



2(n + 5) 



"5 



in 



1 \l (D + 21- 2) 2 + (2r m lC 2 ) z - 8r 2 V m lC 2 
which is Eq. (22) of Ref. [58] if the perturbed mass m\ — and shifting parameter c = 0. 



2r ' 

/ = 0,1,2,3,---), (53) 



(v) We discuss non-relativistic limit of the energy equation (53). When V = S , Eq. 
(23) reduces to a Schrodinger-like equation for the potential 2V(r). In other words, the non- 
relativistic limit is the Schrodinger equation for the potential — 2V e~ r ^ r ° / [l — e~ r//r °] , r = 
a^ 1 . After making the parameter replacements m + Er — > 2m and m — Er — > —E^r in 
Eq. (53) [58], it reduces into the non-relativistic energy equation of Refs. [10,12,13,15,57,59]: 

2 



E 



NR 



a 2 (D + 21 -l)(D + 21 -3)cq 



8m n 



8m a 2 



(2Vo — mi)(2m — m x ) — a 2 (n + 5) z 
(n~+5) 
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S=- 




which is Eq. (23) of Ref. [57] when c and mxare set to zero. It is noted that Eq. (54) 
is identical to Eq. (59) of Ref. [56] for s-wave in ID when the potential is 2V(r), when 
a becomes pure imaginary, i.e., a — > ia and when we set m = 1, mi = and Co = 0. 
Equation (54) can be reduced to the constant mass (m 1 = 0) case in the three-dimensional 
Schrodinger equation: 



which is identical to the expressions given in Refs. [50,52] when the vector potential is 
taken as 2V(r), c = and u = 1 in Ref. [50]. The numerical approximation to the energy 
eigenvalues in Ref. [50] for the last energy equation was found to be more efficient than the 
approximation given in Eq. (19) of Ref. [50]. Taking V = Zae 2 as in [50], we obtain 



For the s-wave (1 = 0), the above energy spectrum is identical to the factorization method 
[23], SUSYQM [12,13] and NU [46] methods. Expanding the energy equation (53) under 



the terms containing the power of (l/m r ) 2 and (r Vo) 4 , we obtain the relativistic energy 
equation 



which is simply Eq. (24) of Ref. [57], where 5 is given in Eq. (54). The first term is the 
non-relativistic energy and third term is the relativistic approximation to energy. 

Now, let us find the wave function y n (s), which is the polynomial solution of 
hypergeometric-type equation. We multiply Eq. (4) by the weight function p(s) so that 
it can be rewritten in self-adjoint form [45,46] 





the weak coupling conditions [(n + 5)/m r ] 2 <C 1 and [V r /(n + 5)] 2 1, retaining only 




(55) 



[u(s)y' n (s)]' + Xp(s)y n (s)=0. 



(56) 



The weight function p(s) that satisfies Eqs. (9) takes the following form 



p(z) = z 



(l-zf, (3 = 25-1 



(57) 
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which gives the Rodrigues relation: 

y nl {z) = B nl z- 2 ^(1 - z )-f>£- [ z n+2 ^{\ - z) n+ ?] 

= B nl P^\l-2z). (58) 

On the other hand, inserting the values of cr(s), 7r(s) and r(s) given in Eqs. (26), (31) and 
(32) into Eq. (10), we get the other part of the wave function 

<f>(s) = z £nl (l-z) 5 . (59) 

Hence, the wave function g n (z) = (f) n (z)y n (z) becomes 

g(z) = C nl z^(l-z) s Pi 2 ^\l-2z) 

= C nl A°\l - z)*P!f™' P \\ - 2z), z G [1,0). (60) 
Finally, the radial wave functions of the Klein-Gordon equation are obtained as 

R x {r) = N nl r^ D -^ 2 e- £ ^ )ar (1 - e~ ar Y P^ ^ (1 - 2e~ ar ), (61) 

with 

(D) _ 1 f ~ , Vc 2 a*{D + 2l-l){D + 2l-3)c ~ 



hca 



moc 2 ) 2 + ^- - El v 



(3 = ^-^4 (S 2 - V 2 ) + (hca) 2 (D + 21 — 2) 2 , 6=±(1 + P), (62) 
where E n i is given in Eq. (37) and N n i is the radial normalization factor. The Jacobi polyno- 

(2e (D \/3) 

mials P n nl ' (1 — 2e ar ) [60] in the last result can be written in terms of the hypergeometric 
function 2 Fi(-n, n + 2e { ^ ) +0+1, 2e% ] ; e~ ar ) which gives the same result obtained in Ref. 
[57]. 

(i) The exact radial wave functions for the s-wave Klein-Gordon equation in ID reduces 
to the following form (in h — c — 1) : 

R n {x) = Cne'V^K* (1 _ e - aS /ro)(l+a)/2p f (2ro V ^r^ ) «) (1 _ ^-x/ro^ 



a = ^Jl + Arl(s 2 -V 2 ), (63) 

where E n is given in Eq. (41). The last formula is identical to Eq. (35) of Ref. [56] when 
the modified scalar potential, So, is se t to zero. 
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(ii) Choosing the atomic units h/2ir = h = c = 1, the exact radial wave functions for the 
s-wave Klein-Gordon equation in 3D reduces to the following form: 



R n (r) = N n e-^^ r (1 - e -Ao)(iH-a)/2 p ^W^-^) (1 _ 2e -r/r 0)? 



pt° Vm ° K ' a \l-2e- r ^) = 2 F 1 {-n,n + 2r oy Jm* - E* + a+l,2r oy Jm? - E^-e~ ar ), (64) 

where _E„ and a are given in Eq. (52) and Eq. (63), respectively. The last formula is identical 
to Eq. (22) of Ref. [57] when the perturbed mass mi is set to zero. 

(iii) The quasi-exact radial wave functions for the /-wave Klein-Gordon equation in 3D 
reduces to the following form (in h — c — 1) : 



Rni{r) = N ni e 



/ m 2 + ^o_ E V (2r Iml+m^tf^) 

V ''o ni (\-e~ r/r °)( +ai)/ P n (l-2e~ r/r °), 

(2r o.k + H^-^a,) 

Pn V o (l-2e" r/r() ) 



= 2 i^-n, n + 2r Wmg + ?&±i^ - E* +a l + l, 2rJml + ^ + - Kv e _ar ), 



a^^(2Z + l) 2 + 4r 2 (^-y 2 ), (65) 

where is given in Eq. (43) and a = r^ 1 . It is identical to Ref. [57] when mi = 0. The 
eigenfunctions in the constant mass case are written as 

Rni(r) = N nl e V (i_ e -rAo)(i+a ; )/2 Pn V '0 (l-2e- r/ro ), 

6| = 7(2i+l) 2 + 4r 2 (Sg-V 2 ). (66) 

At the end of these calculations, the total wave functions of the Klein-Gordon equation with 
position-dependent mass for the scalar and vector Hulthen potentials are 

Vfo^x) = N nl r^ D -^/ 2 e- £ ^ )ar (1 - e~ ar ) s P^ ' P) {I - 2e~ ar ) 



D-2 



eM±ihOi) II ^ 2^+1^+7-2) (sin ^' )l "" J ^-^ +(j - 2)/2 ^-^ + °- 2)/2) (cos^) 
V 2r (n D _x + 2m') {sm ° D - l) P — ( 67 ) 
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where e { nl ' and (3 are given in Eq. (62) and E nl is given in Eq. (37) [46]. 

To check the accuracy of the resulting analytical expressions. We give a few numerical 
real eigenvalues for some selected values of the mass m and mxand potential parameters 
Sq and Vq. In Tables 1 and 2, taking a — 1 and mo = 1, we present some numerical values 
for the energy spectrum of the constant mass Klein-Gordon equations with the condition 
So = Vq for all possible real eigenvalues. To get more real energy eigenvalues in the constant 
mass case (e.g., m = 1, m\ — 0), the vector parameter V of the Hulthen potential should 
be increased. As shown in Tables 1 and 2, when the parameter V — S — 1, 2, 3, 6, 8, 20, 
we obtain N = 1,3,6,10,15,36 real energy eigenvalues, respectively. The numerical solution 
of the position-dependent mass case with vector and scalar Hulthen potential parameters 
satisfying the conditions So = ±Vq and So > Vo are presented in Table 3. For example, in 
Table 3, when the Hulthen potential parameter V — So — 1, m = 5 and m\ ^ 0, we obtain 
N = 46 real energy eigenvalues. Obviously, the number of real eigenvalues increases in the 
solution of the position-dependent case than in the constant mass case where the condition 
So > Vq must be fulfilled. 



IV. COCLUSIONS 



In summary, we have proposed an alternative approximation scheme for the centrifu- 
gal potential similar to the non-relativistic case. This is because the usual approximation 
[10,13,15] for the centrifugal term is only valid for low-screening region, however, for the 
high screening region where a increases, the agreement between the old approximation and 
centrifugal term decreases. Using this approximation scheme, the analytical solutions of the 
radial Klein-Gordon equation with position-dependent mass for scalar and vector Hulthen 
potentials can be approximately obtained for any dimension D and orbital angular mo- 
mentum quantum number I. It is found that the expressions for the eigenvalues and the 
corresponding eigenfunctions become complicated and tedious since the eigenvalues are re- 
lated to the parameters m , mi, So, V , c and a. We have investigated the possibility to 
obtain the bound-state (real) energy spectra with some constraints to be imposed on the 
parameters and, further, the relationship between the strengths of vector V and scalar So 
coupling parameters. In one- and three-dimensions, the special cases for the angular momen- 
tum / = 0, 1 are carried out in detail. We find that the analytical expressions of the energy 
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eigenvalues and eigenfunctions are identical with the results obtained by other methods. 
The analytical energy equation and the unnormalized radial wavefunctions are expressed in 
terms of hypergeometric polynomials. For constant mass case (m 1 = 0) and s-wave (1 = 0), 
the results are reduced to exact solution of bound states of s-wave Klein-Gordon equation 
with scalar and vector Hulthen potentials. To test our results, we have also calculated the 
energy eigenvalues of a particle and antiparticle for the constant mass limit as well as the 
position-dependent mass case. The case of spatial-dependent mass with scalar potential 
So = rriic 2 is found to be equivalent to the constant mass with scalar potential S = in a 
pure vector case. Hence, the spectrum is found to be same. 
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TABLE I: The energy spectrum of the scalar and vector Hulthen potential for rriQ = 1 and m\ = 0. 

V = S nl E+« £-« gj [61,62] fe gjj [61,62] fe 

1 1 1.000000 -0.600000 1.000000 -0.600000 

11- 

2 1 0.707107 -0.707107 0.707107 -0.707107 
1 1 0.984171 -0.214941 - 

1 2 - 

2 0.984171 -0.214941 0.984171 -0.214941 
2 1- 

3 1 0.302169 -0.763708 0.302169 -0.763708 
1 1 0.911438 -0.411438 - 

1 2 0.600000 0.600000 - 

1 3 - 

2 0.911438 -0.411438 0.911438 -0.411438 
2 1 0.600000 0.600000 - 

2 2- 

3 0.600000 0.600000 0.600000 0.600000 
3 1- 

6 1 -0.355051 -0.844949 -0.355051 -0.844949 

1 1 0.235890 -0.635890 - 

1 2 0.763708 -0.302169 - 

1 3 0.994273 0.284416 - 

2 0.235890 -0.635890 0.235890 -0.635890 
2 1 0.763708 -0.302169 - 

2 2 0.994273 -0.284416 - 

2 3- 

3 0.763708 -0.302169 0.763708 -0.302169 
3 1 0.994273 0.284416 - 

3 2- 

4 0.994273 0.284416 0.994273 0.284416 

"The present NU method. 

ft The results from AIM and SUSY. q 



TABLE II: The energy spectrum of the scalar and vector Hulthen potential for ijiq = 1 and m\ = 0. 
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-0.377770 


7 





0.825665 


0.208818 


1 


5 


0.523260 


-0.122370 


7 


1 


0.998229 


0.706553 


1 


6 


0.825665 


0.208818 


8 





0.998229 


0.706553 


1 


7 


0.998229 


0.706553 


9 
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TABLE III: The energy spectrum of the scalar and vector Hulthen potential for mi / 0. 
rriQ m\ Vq So n I E + E~ tuq mi V$ Sq n I E + E~ 



5 0.01 2 2 1 
1 
2 
2 
2 
3 
3 
3 
3 

5 0.01 -2 2 1 
1 

2 
2 
2 
3 
3 
3 
3 

5 0.1 -2 5 1 
1 

2 
2 
2 
3 
3 
3 
3 






0.822925 


-4.913410 5 1 


1 


3.110670 


-4.804170 





3.065630 


-4.807820 


1 


4.252020 


-4.650830 


2 


4.795730 


-4.445800 





4.229630 


-4.655840 


1 


4.793910 


-4.447040 


2 


4.989330 


-4.185200 


3 


4.956220 


-3.857960 





4.913410 


-0.822930 5 0.1 


1 


4.804170 


-3.110670 





4.807820 


-3.065630 


1 


4.650830 


-4.252020 


2 


4.445800 


-4.795730 





4.655840 


-4.229630 


1 


4.447040 


-4.793910 


2 


4.185200 


-4.989330 


3 


3.857960 


-4.956220 





4.871650 


-3.222360 


1 


4.926240 


-3.503700 





5.000000 


-4.245710 


1 


4.995470 


-4.392630 


2 


4.965180 


-4.615030 





4.915250 


-4.768460 


1 


4.878060 


-4.836860 


2 


4.793250 


-4.930300 


3 


4.647670 


-4.993400 



20 1 





4.857570 


-1.483450 


1 


1 


4.875450 


-1.571890 


2 





4.999480 


-2.709050 


2 


1 


4.999990 


-2.772530 


2 


2 


4.998750 


-2.895220 


3 





4.924130 


-3.601650 


3 


1 


4.914310 


-3.648140 


3 


2 


4.893220 


-3.737900 


3 


3 


4.858140 


-3.864780 


1 1 





3.443410 


-4.868720 


1 


1 


4.722690 


-4.742880 


2 





4.618770 


-4.768190 


2 


1 


4.982510 


-4.577550 


2 


2 


4.964780 


-4.347700 


3 





4.960360 


-4.613290 


3 


1 


4.967570 


-4.354450 


3 


2 


4.788530 


-4.056980 


3 


3 


4.484330 


-3.682040 


4 





4.984480 


-4.401670 


4 


1 


4.794830 


-4.065620 


4 


2 


4.488330 


-3.686650 


4 


3 


4.054980 


-3.206920 


4 


4 


3.455290 


-2.575480 


5 





4.837690 


-4.126180 


5 


1 


4.497830 


-3.697630 


5 


2 


4.060510 


-3.212870 


5 


3 


3.459590 


-2.579950 


5 


4 


2.567010 


-1.664550 


5 


5 
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